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Abstract 



We use a generalized 5N formalism to study the generation of the primordial curva- 
ture perturbation in the curvaton brane scenario inspired by stringy compactifications. We 
note that the non-Gaussian features, especially the trispectra, crucially depend on the decay 
mechanism in a general curvaton scenario. Specifically, we study the bispectra and trispectra 
of the curvaton brane model in detail to illustrate the importance of curvaton decay in gen- 
erating nonlinear fluctuations. When the curvaton brane moves non-relativistically during 
inflation, the shape of non-Gaussianity is local, but the corresponding size is different from 
that in the standard curvaton scenario. When the curvaton brane moves relativistically in 
inflationary stage, the shape of non-Gaussianity is of equilateral type. 
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1 Introduction 

As an alternative to inflation, the curvaton model [1] is a well-motivated proposal for explain- 
ing the observed scale- invariant primordial density perturbation in the framework of inflation 
(and some pioneering works on this mechanism have been considered in Refs. [21 El H]). This 
model is based on the inflationary scenario of multiple fields, in which two fields are required 
at least. In this model the universe is composed of radiation decayed from the inflaton 
and the curvaton field. During inflation, the curvaton is subdominant and provides entropy 
(isocurvature) perturbation, and afterwards, the entropy perturbation can be converted to 
curvature perturbation as long as the curvaton decays into radiation before primordial nu- 
cleosynthesis. After the curvaton decays, the universe enters the standard thermal history, 
and then the primordial curvature perturbation leads to the formation of the large scale 
structure of our universe 0, E]. 

Recently, a curvaton scenario realized in the frame of stringy inflationary models was 
presented in Ref. [7J. The idea of brane inflation proposed in jS], is recently successfully 
realized by virtue of warped compactifications[9]. In the stringy landscape a number of probe 
branes are allowed to move in warped throats, and under the description of supergravity 
(SUGRA) their dynamics is described by a Dirac-Born-Infeld (DBI) action in the relativistic 
limit. When applied to cosmology, it was found that a relativistically moving D3-brane in 
a deformed AdS throat is able to drive inflation without slow-roll at early universe |10l [TT] . 
This idea was later extended into the scenario of multiple brane inflation appeared in Refs. 
[T2"j [TBI (see [Tlj [To] for earlier studies under slow- roll approximation) . A distinguish feature 
of the model in this type is that a large positive non-local non-Gaussianity can be obtained 
due to an enhancement of a small sound speed for the perturbation [T6| [TT] . which is in 
contrast to the prediction of a canonical single- field inflation model [18]. Established on the 
scenario of multiple brane inflation [12 [ [T3] and motivated by the model of spinflation [T9] . 
an explicit construction of the curvaton brane model was recently realized in Ref. [20] in 
which angular degrees of freedom had been introduced to perform a relativistic rotation in 
a warped throat. 

The investigation of non-Gaussianities is crucial in distinguishing the curvaton scenario 
from usual inflation models. A key prediction of a general curvaton model is that a sizable 
local non-Gaussianity can be obtained of which the value is mainly decided by the occupation 
of curvaton energy Q x after inflation characterized by a transfer efficiency parameter 
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which is usually calculated on the hypersurface of the curvaton decay. As a consequence, 
the decay mechanism of a curvaton field plays an important role in determining the detailed 
information of primordial perturbation seeded by the curvaton field [5T|. It was pointed 
out in Ref. [21] that at third order, the primordial curvature perturbation under the as- 
sumption of sudden curvaton decay is nicely consistent with that obtained in the case of a 
non-instantaneous decay in the limits of small r and r — > 1; however, their bispectra and 
trispectra deviate from the others in the middle band of the allowed value regime of r. The 
implications of the non-Gaussian features of curvature perturbation on the curvaton scenario 
was studied in [221 [231 [21] . 



In the present paper we extend the investigation of Ref. [21] relating non-Gaussian 
features of primordial curvature perturbation to the curvaton decay mechanism into a more 
general curvaton scenario. Specifically, we consider a model of curvaton brane. We study 
the generation of curvature perturbation and non-Gaussianity originating from the vacuum 
fluctuations on the curvaton brane in cases of a sudden decay on the slice of uniform curvaton 
density, and then compare their observational signatures with the standard curvaton model 
to illustrate the significance of the decay process of the curvaton field. Moreover, since the 
model of curvaton brane involves a sound speed parameter which describes the propagation 
of the curvaton fluctuations during inflation, it has interesting implications to cosmological 
observations as well. Note that in our scenario, the inflaton field can be realized by any string 
realization of inflation which is independent of the curvaton field. Thanks to the virtue of 
curvaton scenario, the detailed realization of inflation does not affect the calculation for 
curvaton perturbation, as long as the inflaton perturbation itself is negligible. 

Our paper is organized as follows. In Section 2, we point out that the non-Gaussian 
features in a curvaton mechanism crucially depend on the curvaton decay mechanism. We 
suggest the curvaton could decay on a slice of uniform curvaton density, and non-Gaussian 
fluctuation evaluated on this slice is mainly determined by a combination of the equation- 
of-state (EoS) of the curvaton and a generalized transfer efficiency parameter which will 
be shown in the main text. In Section 3, we study the model of curvaton brane as a 
specific example. We firstly describe the background evolution of a curvaton brane model 
in relativistic limit throughout the inflationary period and the process of curvaton decay. 
Afterwards, we develop an analytic investigation of the curvature perturbation and the non- 
Gaussianity generated in this model by solving the perturbation equation of the curvaton 
brane order by order. In Section 4 we consider a class of plausible curvaton brane dominant 
eras and derive the nonlinearity parameters characterizing non-Gaussian distribution of the 
primordial curvature perturbation in these cases respectively. Finally, we present concluding 
remarks in Section 5. 

2 Curvaton decay and Non-Gaussianities 

In this section, we first briefly review the standard calculation of bispectrum and trispectrum 
for curvaton mechanism, with the assumption that the decay hypersurface of the curvaton 
field is a uniform total energy density slice. Afterwards, we show that these non-Gaussian 
features behave different under different curvaton decay mechanisms. Especially, when the 
curvaton decay hypersurface is a uniform curvaton energy density surface, the trispectrum 
is significantly modified. We then apply our result to a curvaton model realized by a non- 
relativistic moving brane in a warped throat and study its cosmological implications. 

2.1 Curvaton decay at uniform total energy density 

The non-Gaussian features for curvaton decay at uniform total energy density slice is well 
studied in the literature. Here we briefly review the approach in [2"T] . 

In the local ansatz of non-Gaussianity, the curvature perturbation can be expanded order 



by order as follows, 
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where Ci is the Gaussian fluctuation, and ( n are the non-Gaussian components of order £". 
The relation between ( n and the non-Gaussian parameters yields the following non-Gaussian 
estimators in a general curvaton model, 
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The correlation functions are defined as 
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where P{k\) is related to the dimensionless power spectrum in form of 

Vdh) = £fP{ki) ■ (4) 

To insert the ansatz ([TJ into Eq. (J3J), one can relate the bispectrum P and the trispectrum 
T with P as follows, 

P(k l5 k 2 , k 3 ) = p NL [P{k x )P{k 2 ) + 2 perm.] . (5) 
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+T7VL[P(A;i)P(A; 2 )P(|k 1 + k 3 |) + llperm.] . (6) 

On super-Hubble scales, the curvature perturbation on uniform density slice can be writ- 
ten as 

Qix) = 5N(x) + - / _*..-> , (7) 

3 J Pi {t) Pi + Pi{pi) 

where the subscript i denotes either inflaton, curvaton or total energy density. As there 
is no interaction between inflaton and curvaton, the inflaton perturbation, which will be 
transferred into radiation after inflaton decay, ( r , and the curvaton perturbation ( x are both 
conserved on super-Hubble scales. 

As a first step, one need to write down the relation between the curvaton fluctuation 5x 
and its curvature perturbation. Choosing the spatial flat slice, Eq. ([7]) for curvaton becomes 

Px = Px£ Kx , (8) 



in the neighborhood of the curvaton decay hypersurface. 

Consider the curvaton perturbation that is initially originated from vacuum fluctuations 
inside the Hubble horizon. These perturbation modes could satisfy a Gaussian distribution 
at the Hubble exit. For a curvaton brane, this assumption applies when the brane is moving 
non-relativistically. In this case, we have 

X* = X* + ^iX* • (9) 

This Hubble-crossing value can be related to the initial amplitude of curvaton oscillation, 
which takes the form 

t^i nl V 9' J 

The detailed form of g(x*) is a mo del- dependent function which is determined by the explicit 
potential of the curvaton field in the non-relativistic limit. For example, if the curvaton 
potential is quadratic all the way until curvaton decay, then g(x*) oc x*- I n this case, the 
energy density can be written as 



2 2 ^2 

p v = —m q = -m 
Hx 2 2 



+ yy_rnx 

n=l 

Therefore, the combination of Eqs. ([H]) and ( TTT1) leads to 
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A second step of the calculation is to relate ( x to (. In the sudden decay approximation, 
the relation is quite simple and can be computed analytically. However, we also note that 
this relation crucially depends on the curvaton decay mechanism, which will be discussed 
in the next subsection. In this section, we assume the curvaton decays on a uniform total 
density hypersurface H = T, where T is the decay rate of the curvaton. This can be realized 
in the framework of particle physics, namely the curvaton has a "life time" and then decays 
to particles through their coupling terms. Then on the curvaton decay hypersurface we have 

Pr + P x = P ■ (13) 

In a usual case the curvaton is oscillating around its vacuum before its decay. Thus one 
obtains a pressureless field fluid on average with P x = for curvaton, and P r = p r /3 for the 
radiation from inflaton. Making use of Eq. ([7|), one can get 

Pr = p r6 ^-0 j px = p xe ^-0 . (14) 

As a consequence, ( and £ x are related on the decay hypersurface as follows 

(1 - l\)e 4 ^-« + n x e 3 ^~V = 1 , (15) 



where fi x = p x /(p r + P x ) * s the dimensionless density parameter for the curvaton at the 
decay moment. A necessary condition for the curvaton mechanism working is based on the 
assumption that the fluctuation ( r seeded by the inflaton field is negligible. Then ( and C x 
can be related order by order as 

Q. 2r\ + g 
Cs 9 g(3g'g" + gg'") 9 / gg"\ 9gg" 1 1fl ^„ 2 ,„. 

where r is the transfer efficiency defined in the Introduction. Finally, we have the nonlinearity 
parameter Jnl in a usual curvaton model as follows, 
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Additionally, the non-Gaussian estimator of the third order fluctuation can be simplified in 
the absence of nonlinear evolution of the curvaton field between the Hubble exit and the 
start of curvaton oscillation which yields g" = g'" = 0, 
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Notice that Eq. (!T8|) corresponds to the simplest example where a canonical curvaton has a 
quadratic potential. In this case, one obtains g^L oc r ~ 1 which indicates that the trispectrum 
is more difficult to be detected in the simplest curvaton model. However, if g is of a nontrivial 
form, the trispectrum g^i could scale as r~ 2 in the small r limit. 

We would like to emphasize two key features of the above calculation. Firstly, the curva- 
ture perturbation (, ( x and ( r studied in the above are gauge invariant. Therefore, although 
we study their dynamics in the gauge of spatial flat slice, their values are irrelevant to the 
gauge choice. Secondly, in the sudden decay approximation, ( is conserved right after the 
universe evolves through the curvaton decay hypersurface and the curvature perturbation 
should be calculated exactly on this hypersurface. This is because, if one calculates at any 
earlier time, ( is not conserved; if one calculates at any later time, one can no longer use the 
information stored in ( x . Having the above comments in mind, it is clear that if the curvaton 
decay surface is different from the uniform total energy density surface, the non-Gaussian 
features are also changed. This is the right topic we shall investigate in the next subsection. 

2.2 Curvaton decay at uniform curvaton density 

In this subsection, we phenomenologically assume that curvaton decay happens on a uniform 
curvaton density slice. The realization of this assumption on a brane is discussed in the next 
subsection. 

If curvaton decay does not happen on the uniform total energy density slice p = p(tdec) 5 
instead, happens on the uniform curvaton energy density slice 

Px = Pxitdec) , (19) 



we can calculate the relation between £, £ r and C x on ^he uniform curvaton density slicqj. 
In this case we apply Eq. ([7j) on the slice (I19p and thus have 

( X = 5N . (20) 



For radiation, we have 

t. = C. + 

and from Eq. (1211) . we can solve p r as 



C, = C x + 7 lo g(?) > (21) 
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Pr = p r e 4 ttr-<x) . (22) 

For the total energy density, we have 



^ = ^ +i °g ; xp ' (23) 

4 VPr + 3(p x + P x )/4y 

in which we have applied the relation dp = dp r , since on the slice of uniform curvaton density 
p x is a constant and can be absorbed into the background. 
Inserting Eq. ( 12 2 p into Eq. (I23p . we have 



1 _ ±_^n, e 4 «"« = (1 - a )e 4 «-"« + k I ; »v , (24) 



where we have introduced an EoS parameter u> = P x jp x for the curvaton. This equation 
should be used to replace Eq. (115]) for doing the subsequent calculation when the curvaton 
decay slice has uniform curvaton energy density. 

Again, we consider the limit of negligible ( r , as is assumed in a usual curvaton scenario. Of 
course, one could generalize the calculation into the case of mixed perturbations of infiaton 
and curvaton, however, this is beyond the scope of the current paper. Thus the relation 
between ( x and 5\* is in principle not modified compared with the last subsection. However, 
if we want to consider general EoS for curvaton, the relation in the spatial flat slice (analog 
to Eq. (JSD) becomes 

Px = P x e 3il+w)Cx ■ (25) 

When 5x* is Gaussian, to expand the above equations order by order, we have the following 
generic relation for linear curvature perturbation: 

Ci = rC x i , (26) 



Another approach to calculating curvature perturbation can be chosen on the slice of uniform total 
density at some moment 4y after curvaton decay. At that moment, we have the total energy density ptotitr) = 
Pr, x {t T )e i{c ^~ SNT) + p r (t T )e^ Cr T-^ T ) anc } the curva t U re perturbation ( T = SN T , where the subscript "T" 
denotes the time tx- However, after curvaton decay the universe has already entered a radiation dominant 
stage and thus is always adiabatic which yields a conserved curvature perturbation £x = C{tdec) at large 
scales. Therefore, the result will not change if we compute the generation of primordial perturbation at the 
slice of curvaton decay as shown in the main text. We are grateful to Misao Sasaki for pointing out this 
approach. 



where we have introduced a generalized transfer efficiency as 

(27) 



3(i + w)n x 
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which encodes the information of the EoS of the curvaton field. 

To proceed, we solve the curvature perturbation order by order again and derive the 
non-Gaussian fluctuations as follows, 

|^( 1 + f Q + I^_ 4 , (28) 

Cs 9(1 + wf g(3g'g" + gg'") 9(1 + w)(l - 3w) gg" (5-3w)(l-3w) 
Cf 4f 2 g' 3 2f 2 g' 2 2P 

_18(l + ro )/ 1 + ^_12(l-3 ro)+32 
r \ g' 2 J r 

Substituting the above results into the non-Gaussian estimators (j2j), we are able to obtain 
the nonlinearity parameters of local shape. 

Very interestingly, large trispectrum can be generated even in the case of a constant g. 
This conclusion is quite different from that obtained in the previous subsection. Another 
thing worthy to note is the generalized transfer efficiency parameter f contains a factor 
1 + w in its expression. As a consequence, when w — > — 1, the value of f can be strongly 
suppressed and thus the amplitude of non-Gaussian fluctuations can also be enhanced. This 
result indicates that the non-Gaussianity could be amplified if there exists a secondary 
inflation. This property was earlier discovered in several specific curvaton models [25] and 
[13], in the case that curvaton decay still occurs on the slice of uniform total energy density. 
For example, in [25] a secondary inflation was achieved in the limit of n — > in their model; 
while in [13] the similar background solution was obtained due to the survival of a brane 
with a light mass term. In our work, however, we point out clearly that the mechanism of 
enhancing the non-Gaussian fluctuations due to a secondary inflation is a generic feature in 
a curvaton model if the field decays on the slice of uniform curvaton density as we analyzed. 

2.3 Application to models of non-relativist ic curvaton branes 

The decay mechanism in the previous subsection has interesting application in the case of 
brane inflation. This is because a single D-brane is a BPS object, which is stabilized by 
supersymmetry. The decay of a D-brane is typically expected to be through the collision 
process with an anti D-brane. 

We first consider the case when the curvaton is the position modulus of a D3 brancl 
The annihilation of a curvaton D3 and an anti D3 takes place when the distance between 



One should be careful in imposing initial conditions to a curvaton brane. Usually one can impose 
initial condition in two ways in curvaton scenario. One possibility is that the curvaton originally sits at the 
minimum of the potential. If the total inflationary e-folding number is much larger than 60, the random 
walk of curvaton will finally set the initial condition for curvaton for the last 60 inflationary c-folds. Another 
possibility is to impose the initial condition without the above random walk assumption. For brane inflation, 
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the curvaton D3 and the anti-D3 is of order string length. After that, the potential of the 
brane-anti-brane system becomes tachyonic, and the system decays into lower dimensional 
solitons. 

In the framework of moduli stabilization, it is natural to assume that the position modulus 
of the D3 is stabilized on the tip of the deformed warped throat. This indicates that the 
open string modes characterizing the position of the anti-D3 are very massive, and so the 
anti-D3 is not allowed to move and even its fluctuations are strongly depressed. In this case, 
the relative position between the curvaton D3 and the anti-D3 is described by the position 
of the curvaton D3, up to a constant. Thus the curvaton decay surface becomes the uniform 
curvaton energy density slice, instead of the uniform total energy density slice. 

When the curvaton D3 moves non-relativistically, the brane position can be described by 
a scalar field with a canonical kinetic term. The most natural curvaton decay mechanism 
is that the D3 moves towards the anti-D3 and then annihilates into a system of tachyon 
condensate [32] without an oscillation period. In this case, the EoS of the position modulus 
evolves from w ~ — 1 into w ~ gradually as happened in tachyon cosmology [33]. Thus 
from Eqs. (128]) and (T29]) . the non-Gaussianities can be enhanced during the occurrence of a 
secondary tachyonic inflation with w ~ — 1 which will dilute other matter components in a 
few e-folds; otherwise, we can also obtain large non-Gaussianities due to a small value of Q x 
if the curvaton brane decays into a pressureless tachyon condensate directly. 

Now we consider another curvaton brane model achieved by a probe anti-D3 with angular 
motion at the bottom of a warped throat with approximate isometries [20]. As usual, if a KS 
throat is isometric, its warp factor is independent of the angular coordinates. However, in a 
general case, there are some corrections which make the warp factor dependent on angular 
coordinates. For example, since the compact C-Y manifold cannot have exact continuous 
isometries, the isometries of the bulk must be broken when the finite throat is glued on 
this bulk, and then the warp factor dependents on angles. Moreover, a nonperturbative 
effect which stabilizes the Kahler moduli could bring a potential for the branesJSH [35J, [36] . 
The form of this potential depends on the precise embedding of the wrapped branes. For 
a general embedding which does not admit supersymmetric vacua on the tip, this potential 
could become dependent on angles as well. Namely, an estimate of the nonperturbative effect 
which is dominated by the warp factor in the case of the Kuperstein embedding|37j of the 
D7-brane gives rise to a mass term for the angular degrees of freedom, which is determined 
by the minimal radial location reached by the D7 and the deformation parameter of the 
conifold. 

In the model of curvaton brane involving angular degrees of freedom, one may choose the 
flattest angular direction along the rotation of the curvaton field. Since in this model the 
curvaton brane is acted by an anti-D3, its radial coordinate is almost fixed on the tip of the 
throat as explained above. In the limit of a non-relativistic anti-D3, one could finely tune 
the model parameters to allow the curvaton anti-D3 to rotate slowly so that it can survive 
during inflation. After inflaton decay, the anti-D3 is able to oscillate due to the existence of 



one can show that the quantum fluctuation during one Hubble time is typically smaller than the classical 
motion |26j . If this conclusion also applies for curvaton, one should expect only the latter possibility for 
brane curvaton initial condition. On the other hand, our scenario is motivated by string theory. As inflation 
in string theory usually needs some fine-tuning [27], one expects the e- folding number is not much larger 
than 60. This also supports the latter possibility. We thank Xingang Chen for discussion on this point. 



a nearly quadratic potential in angular space. 

In the above, we have introduced two explicit realizations of the curvaton field in string 
compactification, i.e., by the distance of a light D3 brane in the throat, as well by the flattest 
angular degree of freedom of anti-D3 stabilized at the bottom of the warped throat. For a 
light D3 brane, it decays only when it arrives at the tip of the warped throat and annihilates 
with an anti-D3; while for the angular degree of freedom of an anti-D3, it decays only when 
the anti-D3 rotates to the position which is parallel to others. In these cases, the value of the 
curvaton field is almost the only parameter which determines the annihilation of curvaton 
brane, and thus the curvaton decay. In this sense, the inflaton energy density is not relevant 
to the curvaton decay. Therefore, the curvaton decay slice is the uniform curvaton field value 
slice, and thus the uniform curvaton energy density slice. 

Finally, we would like to point out that our calculation of curvature perturbation on uni- 
form curvaton density decay hyper-surface is not limited to curvaton brane, but is generically 
applicable to any other models in which the curvaton decays in the uniform curvaton density 
slice. For example, our calculation could be applicable for the curvaton scenarios in which 
the curvaton decays from a narrow band preheating. In this preheating type of curvaton 
scenario, curvaton decays each time when its field value crosses zero, which corresponds to 
a uniform curvaton energy density slice instead of a uniform total energy slice. Our method 
also has significant implications if applied to the curvaton scenario in the frame of bouncing 
cosmologjo However, since these scenarios are beyond the scope of the current work, we 
would like to leave the detailed analysis on these issues in future projects. 

3 The model of a relativistic curvaton brane 

In this section we extend the approach to studying the generation of non-linearities on the 
uniform curvaton density slice into the case involving a sound speed parameter. This can 
be obtained in an example of a relativistic curvaton brane model. Specifically, in the frame 
of warped compactifications, we consider a system constructed by a number of D3-branes 
in a background metric field with negligible covariant derivatives of the field strengths and 
assume that these branes are falling into the warped throats, this system can be described 
by a DBI action which has a stringy origin. 

3.1 Background analysis 

We phenomenologically consider a double field inflation model with its action in form of 

C = V=9 U - V(<f>) -ly/l-2fY- W( X ) } , (30) 

where 



5 Some observational signatures of bouncing cosmology were extensively studied in the literature, such as 
curvature perturbations [28l 129) . non-Gaussianities|30|. and constraints from current observations |31j. 



We shall find a solution in which behaves as inflaton, which dominates the energy density 
during inflation, and \ behaves as curvaton, which can be realized by a mobile probe brane 
in a warped throat during inflation. 

To assume the curvaton \ is moving relativistically, we can obtain a small sound speed 
defined by 



c s = VI - 2/y < 1 , (32) 

at the limit of 2fY — > 1. The continuous equation of x leads to 

_y± — j_x_ w , (33) 

A 3H K J 

One can solve \ 2 from this equation as 

2F = v 2 = ^ 3g ) (34) 

To make this equation consistent with the inequality ( 13"2"j) . the following condition is expected 
to be satisfied in the model, 

fW'/(3H) > 1 . (35) 

Another necessary condition is that the assumption that inflaton dominates the energy den- 
sity: 

V > W . (36) 

Making use of Eq. ( 13 6 p and the Friedmann equation 3M 2 H 2 ~ V, Eq. ( 13 5 p can be rewritten 
as 

f 2 W' 2 M 2 > 3V > 3W . (37) 

Assuming that the curvaton potential is a mass term W = m 2 x 2 /2, we need 

fm 2 Ml » 1 . (38) 

If this condition is satisfied, we could have certain parameter space to let x behave as a 
relativistically moving curvaton. Theoretically, this relation could be possible in condition 
of that the warp factor h is finely tuned to be sufficiently small by adjusting the position of 
the curvaton brane in the throat. In this case a sufficiently small warp factor yields a large 
value of / and therefore the field x is a ble to slowly roll down along the potential even this 
potential is not enough flat. Consequently, x is a ble to stay at sub-Planckian regime safely, 
so that the backreaction on the warped throat could be controllable. 

Now let us examine the dynamics after inflation. Depending on the choice of / and m, 
as well as the brane position, the curvaton might keep to stay in the relativistic regime, or 
slow down and approach to a standard kinetic term. Practically, if the curvaton keeps in the 
relativistic regime until it decays, there is no oscillation stage of curvaton, and a secondary 
inflation driven by curvaton is possible when the curvaton start to dominate. On the other 
hand, when the curvaton becomes non-relativistic before its decay, one can either have a 
stage of secondary inflation or curvaton oscillation. 
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3.2 Perturbation 

In a relativistic version of curvaton brane scenario, the calculation of perturbation is the same 
as that of a single scalar field with generalized kinetic terms in the slow roll approximation. 
The only change is to write ( x as the perturbation variable instead of (. In this section, 
we shall first calculate the three point function explicitly, which shows the equivalence with 
single field DBI inflation. After that, we directly apply the result of four point function of 
DBI inflation to a relativistic curvaton brane scenario. 

In the curvaton scenario, one can neglect the perturbation of the inflaton field. As the 
curvaton field is subdominant during inflation, one can ignore the perturbation of the metric. 
The second order action of the curvaton takes the form 

in which the next-to-leading order terms have been neglected due to a suppression of the 
slow roll condition. 

We use C x to denote the curvature perturbation in the uniform curvaton energy density 
slice. To leading order in c s , one can write 

Cx = -^- (40) 

X 

We assume that the curvaton field is not directly coupled to the inflaton field (except for 
the coupling via gravity). In this case, as shown in [3H], Cx * s a conserved quantity on super- 
Hubble scales. At the moment of Hubble-crossing, one can solve out 5x* = IpS and so the 
corresponding power spectrum is given by 

P C = ? 2 Pt x = r^ 2 . (41) 

3.2.1 Bispectrum 

To expand the action into cubic order, we obtain the dominant part of the action as follows, 

3 

2 V c s J \ c3 s X a 2 c s x 

The initial correlation function of x can be calculated using the time dependent perturbation 
theory. Firstly, the mode function of the perturbation can be written as 

u k = -jL= (1 + ikc s r) e~ ikCsT , u' k = I ^c 2 s re- ikc ° T . (43) 

One can show that the first term in the third order action leads to a contribution to the 
three-point function as 

„ Jff" / ' A M'f'Ek.) 



SlWJ i4(' ^a.w, (42) 



«*> - w U " V ^mMr • (44) 



ll 



and the second term in the third order action leads to a contribution to the three-point 
function as 

. 3 , _ H» (\ \ (27r)W(Ek,)(k 2 .k 3 ) / 1 k 2 + k 3 2fefe 



(c *> 3 ~w U " V %w — L U + -*r + ~w) + 2 penn - • (45) 

Summing them up, we have 



Substituting Eq. ( 1261) into Eq. (1461) . the three-point function for curvature perturbation (, 
takes the form 



(C 3 ) 

r 



? (? - p l ( 2 *) T * 3 E k ') f -jr E *. 2 *f + 5^ E *W + 1 E *f ) + «v 



(47) 



where K 3 = ki + k 2 + k 3 , and (C 3 )ioc is given by Eqs. (J5J) and (J2J). The leading order non- 
Gaussianity is a combination of equilateral shape and local shape. In the c s — > limit, the 
shape becomes completely equilateral. When c s is not very close to zero, the contribution 
(C 3 )ioc is also important. 

3.2.2 Trispectrum 

For the trispectrum, as shown in [39j HO], the correlation function includes a scalar propa- 
gation part and a contact interaction part. The detailed shape is complicated. Instead of 
providing the complete form of the lengthy shape function (as shown in the appendix), we 
would like to study the representative shape here. The trispectrum can be approximated by 

(C 4 ) = 4( 2 -) 9 ^^ 3 (E k ^4^P + KV + (Ocross , (48) 

r c s {i\a) 

where if 4 = k\ + k 2 + k 3 + fc 4 , and a is a constant in order of unity. The local part (C 4 )ioc is 
given by Eqs. (jHj) and (EJ. The cross correlation can be calculated as 

(C 4 )cross =§/&<(C * C)(ki)C(k 2 )C(k 3 )C(k 4 )) + 3 perm. 
5 

= l/j^PCfeXCCki - k 2 )C(k 3 )C(k 4 )> + 11 perm. , (49) 



and (C(ki — k 2 )C(k 3 )C(k 4 )) is calculated using Eq. (I4"6j) . 

4 Non-Gaussian features of curvaton fluctuations in 
specific examples 

In the above section we have studied the generation of curvature perturbation up to third 
order. In order to obtain a much intuitive insight of these non-Gaussian features, one still 
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needs to evaluate the nonlinearity parameters f^L and g^i in explicit cases. As we have 
introduced in Section 2.3, there exist two plausible evolution trajectories for the curvaton 
brane after inflaton decay. We study these possibilities respectively in the following. In 
these specific examples, we generically consider the curvaton potential is almost quadratic 
around the moment of curvaton decay so that we are able to apply the formalism developed 
in Section 2. Therefore, when we compute the non-Gaussian fluctuations at curvaton decay, 
the probe brane has already become non-relativistic as a canonical field. However, one may 
notice that the sound speed parameter c s still exists in the following computation, since 
it was inherited from the Hubble-crossing of curvaton fluctuation during inflation. On the 
other hand, when the probe brane is still moving relativistically near curvaton decay, one 
can use the techniques developed in [13] to generalize the calculation of this Section. 

4.1 Case of w — > — 1 

We first consider the case that near the curvaton decay, the curvaton brane has EoS w — > — 1. 
The theoretical realization of this scenario was discussed in Ref. [13]. If we consider the 
curvaton brane is falling into the AdS-like throat and neglect the backreaction of the branes 
upon the background geometry, the warp factor usually takes the form of f(x) — -Vx 4 - 
Note that, the amplitude of the curvature perturbation depends on the term \ of which the 
lower bound corresponds to its value near the bottom of the warp throat. The warp factor 

at the infrared (IR) end of the throat is given by hm ~ exp ( —j 11 ^ ) where g s is the string 
coupling, M and K are the RR charge on the S% cycle and NS-NS charge on the dual cycle 
respectively. Therefore, at the relativistic limit, one obtains 



Xm-X 2 m\/^--VT3h 2 IR , (50) 

at the IR cutoff. 

The key parameter related to the generation of primordial curvature perturbation in 
the curvaton model is the perturbation transfer efficiency f defined in Eq. ( I27p . From the 
definition, one learns that this parameter is determined by the EoS w and density parameter 
Q x of the curvaton. If the curvaton brane decays when its energy density is still sub-dominate, 
then there could be a double suppression on the value of f because of a combined effect of 
w — > — 1 and Q x — > 0. 

However, if the curvaton does not decay until it become dominate, the universe enters a 
secondary inflation driven by the curvaton. The energy density of the curvaton will catch up 
with that of radiation decayed from inflaton in the first several e-folds. As a consequence, it 
yields Q x — > 1 at the moment of curvaton decay. Note that Eq. (|27|) can be rewritten as 

(51) 



4(1 - o x )/(i + w) + n x ~ 4(1 - n x )/(i + iu) + i ' 

and thus one concludes that the order of magnitude of f is determined by the ratio between 
1 — Q x and 1 + w. Namely, when (1 — Q x )/(1 + w) 3> 1, one has f < 1; and when 
(1 — Q x )/(1 +to)<l, one has f ~ 3. 
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The correlation functions are given by Eqs. ( 147|) and (148|) respectively. One can decom- 
pose the non-Gaussian estimator into local and equilateral ones by comparing their squeezed 
limit. The local components for f NL and g NL are given by Eqs. (I2"gj) . (12"9"j) and ([2]). The 
equilateral components for j'nl and qnl are given by [T71 [39] 

Jnl - 108f Vc2 V ' 9nl ~ <%* ' ^ ] 

which implies the equilateral non-Gaussian fluctuations can be doubly enhanced if both f 
and c 2 have small values. 

4.2 Case of w — >• 

After inflation, the curvaton brane will arrive at the tip of the throat. At that moment, 
either the open string modes of the curvaton brane lead to a matter-like tachyon condensate 
with w ~ — c 2 — > as studied in Ref. [7], or a non-perturbatively moduli stabilization 
brings a quadratic potential for the angular modes of the brane which gives rise to a non- 
relativistic oscillation of the curvaton with w ~ in average and c s ~ l[41j. Consequently, 
we have two plausible ending processes of the curvaton brane which yields different results 
for non-Gaussianities. 

4.2.1 Local Non-Gaussianity 

According to the analysis performed in Section I2T21 we learn that the non-Gaussian fluctua- 
tions in local limit are insensitive to the sound speed c s but mainly depend on the parameter 
of transfer efficiency f. When w ~ the form of f coincides with the usual transfer effi- 
ciency r. Therefore, following Eqs. (|28|) and (|29|) . the non-Gaussian estimators ([2]) give the 
nonlinearity parameter of local shape as follows, 



5_ L + gf\ + 5_ _ 10 

4r V g' 2 J 6r 3 



/ JV£ ~_(l + ^-)+_-_, (53) 



when the curvaton decay takes place on the slice of uniform curvaton density. 

In the limit of r — > 1 when the curvaton dominates the universe before it decays, the 
nonlinearity parameter is given by, 



5 (l - 9j" 
4 V ff 



/^-►-7(i-*Sr) . (54) 



of which the value is negative and sizable if g" is negligible. On the other hand, we have a 
local non-Gaussianity in the limit of r <^ 1 which takes the form, 



4r \3 g 



fNL^i~A- + y -^) , (55) 



which is amplified by a small r. We notice that the above result is consistent with the result 
obtained in a usual canonical curvaton model only in the case of r — > 1, but deviates from 
the usual model in the case of r — > with a factor of order 0(1). 
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In the absence of higher nonlinear evolution of the \ field between the Hubble exit and 
the decay moment, we would have g'" = for simplicity. In this case, the third order 
perturbation gives 

25(5 -8r)gg" 125 125 400 

9nl * 24r2 ~g^ + 108H ~ ^ + YT ' (56) 



Again we consider two limited cases and list the results as follows, 

(57) 



_ 25 3gg" 

9nl ->- t^(! -TT^" 



-> 



12 v 2 ff 

125 / 9gg" 



108r 2 V 2g 



a 



(51 

r-+0 



The result in the first limit is in agreement with that of a usual curvaton model. However, 
the latter one is quite different since its form is enhanced by order of r~ 2 even when g" = 0, 
but in a canonical curvaton model with a quadratic potential we have gNL oc r _1 when g" is 
negligible. 

In order to understand the above analytic study more clearly, we compare the nonlinearity 
parameters /tvl and gNL of local shape in the model of curvaton brane with the results 
obtained in a usual scenario in Figure. [TJ In this figure, we numerically plot the nonlinearity 
parameters as functions of the transfer efficiency r, and consider the simplest curvaton model 
without nonlinear evolution of the \ field between the Hubble exit and the moment of 
curvaton decay which gives a linear function of g(x). The red solid curve represents for the 
non-Gaussianity of local shape when the curvaton field decays on the hypersurface of uniform 
curvaton density; while, the blue dash line denotes the non-Gaussianity of local shape when 
the curvaton decay takes place on the slice of uniform total energy density. One notices that 
these curves coincide at the limit of r — > 1, but possess different behavior when r is less than 
unity. The result for g NL is more manifest than for f NL . 

4.2.2 Equilateral non-Gaussianity 

Now we study the non-Gaussian fluctuation of equilateral shape in the case of to — > 0. 
Since we have pointed out that the sound speed c 2 s can be either very small in the case of 
tachyon condensate or approach unity when the curvaton obtains a quadratic potential, we 
phenomenologically treat the sound speed as a free parameter in the following. Again, we 
obtain the equilateral components for f NL and g^L in form of 

Jnl = ino ( ~j — 1 J > 9nl — A 9 ' ' ) 

108r \c; J & s r l 

This result also indicates that both small r and c 2 s are able to amplify the equilateral non- 
Gaussianity in the model of curvaton brane with w ~ 0. 

From above analysis, if the curvaton brane moves non-relativistically during inflation with 



c 



,2 



1, there are f^L oc 1/r and g^L oc 1/r 2 for both the local and equilateral shapes, but 
the amplitudes of local shape is much larger than equilateral ones; however, if the curvaton 
brane moves relativistically with c 2 s ~ in inflationary stage, the amplitudes of f^L and 
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Figure 1: Absolute values of the nonlinearity parameters /tvl and g^L of local shape as a 
function of the transfer efficiency r in the model of curvaton brane and the usual curvaton 
scenario. Here we assume that g(x) is linear for simplicity. 

g^L could be greatly amplified by the order of l/c 2 s and \jc\ respectively. As a consequence, 
one concludes that when the mobile curvaton brane is non-relativistic, the shape of non- 
Gaussianity is mainly of local type with its corresponding size being different from that in 
the canonical curvaton model, but becomes of equilateral type in the case of a relativistic 
curvaton brane. 



5 Conclusion 

In this paper, we have calculated the nonlinear primordial curvature perturbation in the 
curvaton scenario using a generalized SN formalism. If curvaton decay does not occur on 
the slice of uniform total energy density but on the slice of uniform curvaton density instead, 
we find that the dynamics of the nonlinear fluctuations behaves different from the usual 
scenario. Specifically, we consider a model of curvaton brane which can provide a theoretical 
realization of the curvaton decay we expected. In the frame of this model, we have presented 
a full analysis on the sizes and shapes of its bispectra and trispectra respectively. We arrive 
at an important conclusion that the generation of nonlinearities is sensitive to the mechanism 
of curvaton decay, which is mainly determined by the EoS w, the density occupation Q x of 
the field x & t curvaton decay and its sound speed c s during inflation. Explicitly speaking, our 
results show that the nonlinearity parameter characterizing the second order perturbation 
Jnl is proportional to 1/f in which the transfer efficiency relies on the EoS w and the density 
occupation Q x of the curvaton. For the equilateral shape, the size of Jnl is amplified by a 
large value of the factor 1/c^, and therefore when the curvaton brane moves relativistically 
the shape of non-Gaussianity is mainly of equilateral type; however, at the squeeze limit 
the nonlinearity parameter f^L decouples from the sound speed and its size is qualitatively 
consistent with the result obtained in a usual canonical curvaton model. Moreover, for the 
third order fluctuations, we find the nonlinearity parameter g^i is enhanced by the factor 
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1/Cjj in equilateral shape and decouples in local limit, and in the simplest example its size is 
proportional to 1/f 2 which is quite different from the usual scenario. 

As an end, we would like to highlight the importance of our study in this paper. The 
mechanism of curvaton decay after inflation could determine the relation of the primordial 
curvature perturbation and the decay hypersurface. This process is rather robust and should 
be considered in any specific curvaton model, and the example of curvaton brane analyzed 
in the current work is a good illustration to emphasize its importance. Additionally, the 
signatures of curvaton decay imprinted on the primordial curvature perturbation could pro- 
vide a new window to explore the combination of early universe physics and astronomical 
observations. 
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Appendix: Full shape function of the trispectrum 

The shape of trispectrum for a relativistic brane curvaton is the same as that when the 
brane behaves as an inflaton. For completeness, here we present the full shape function of 
the trispectrum, which is calculated in [39]. As usual, we define 

2 P,x 

c - 



P x + 2XP r xx 
E = XP X + 2X 2 P XX , 

2 

\ = X P y xX + 77 AT P^XXX , 

H = \x 2 P,xx + 2X 3 P X xx + \x^P,xxxx • (60) 

There are two classes of contributions to the trispectrum, namely, the contact interaction 
and scalar exchange diagrams. They are shown in Fig. El^a) and 121(b), respectively. Using 
the method of time dependent perturbation theory, the former can be calculated from the 
forth order action, and the latter can be calculated by applying the third order action twice. 

We will skip the lengthy calculation, and refer the readers to [39] for details. It turns out 
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Figure 2: Two classes of contributions to the trispectrum. The contribution (a) is denoted 
by "contact interaction" diagram, and the contribution (b) is denoted by "scalar exchange" 
diagram. 



that the trispectrum for (^ can be written as 



4 4 

(0 = (2n) 9 PlS 3 (Y,^)\lT3r(k u k 2 ,k 3 ,h,k 12 ,k u 



i=l i=l 



where T has six components: 

2 



r 



fO' T * + *{h.- 1 ) T * 



1 T, 



s3 



li 9A S 

E ~ E 2 



T, 



c\ 



(61) 



(62) 



where T cl , T c2 and T c3 are contributions from contact interaction, which are given by 

> 9A 2 \nlA 2 



T^36( § --l4^, 



(63) 



1 /3A 



-,:■! 



8 V E 



fc?fc|(k3 • k4) 

K 3 



1 3(fc 3 + fet) I2fc 3 /C4 



A' 



A 2 



+ 23 perm. 



(64) 



- c3 



iri-i 



32 

+ 23 perm. . 



(k 1 -k 2 )(k 3 -k 4 ) 
A 



E 



:j kikj 3fcifc 2 fc 3 fc 4 y^ 1_, 



A 2 



A 3 



i=l 



12 



hk 2 k 3 k 4 



A 4 



(65) 



Tsi) 7s2 an d T s3 come from interaction with scalar exchange. As shown in Eq. (1421) . there 
are two dominant terms in the curvaton third order action. One can show that the coupling 
coefficients of these two terms are of order A/E and \ — 1, respectively. T s \ denotes the 
contribution that both interaction vertices in Fig. EJ^b) have coupling A/E, and thus is 
proportional to (A/E) 2 . T s2 denotes the contribution which is proportional to A/E(4j — 1). 

Additionally, T s3 denotes the contribution with two powers of \ — 1. 



The summation T sl + T s2 + T s3 are given by a summation of the following contributions 
from Eq. (1S"6"1) to Eq. 



| g) ' ^^^ (fei + fc2 ; M 3 M3 + 23 perm. . (66) 



I (|) ' tfM^fc, jj)i (^ + ^ + 7^)+ 2 3 P— • (67) 



4^ ( \ - A < k 3 ' U)k^ikl *_,, , 3 F(k 3 , k t , M) + 23 perm. . (68) 

32L\c| / (A:i + fe2 + fei2) 



^ G? ~ (kl2 ' k4)M ^ WT^Tk^ F ^ **■ M) + 23 perm - • (69) 



^(^-i 



' 1 (k 12 ■ k 4 )^% Gab{k ^ h) + 23 perm> _ (n) 



S \C? / fci2 

~ (^ - l) (ki • k 2 )fc 3 2 A; 2 A; 12 G^fa, fc 2 ) + 23 perm. . (72) 



3 A / 1 \ PPP 

££ (4 - X H k * ■ ki^)^ 1 G ba (-ki2, h) + 23 perm . (73) 

8Y. \cj J k 12 

1 f 1 - 1 J (kx • k 2 )(k 3 • k 4 )fc 12 F(fci, fc 2 , fc t + k 2 + k 12 )F{k 3 , k A , M) + 23 perm. . (74) 

1 / 1 \ 2 k 2 

— ( - - 1 J (ki • k 2 )(k 12 ■ k 4 )^- F(k u k 2 , h + k 2 + k l2 )F{k 12 , h, M) + 23 perm. . (75) 

2 1,21,2 



1/1 \ FF 



- ^ - - 1 j (k 12 • k 2 )(ki2 • k 4 )-^ F(fc 12 , k 2 , h + k 2 + k 12 )F(k 12 , h, M) + 23 perm. .(76) 



1 (\ ^ 2 



— (— -1) (ki -k 2 )(k 3 -k A )k 12 G bb (ki,k 2 ,k 3 ,k 4 ) + 23 perm. . (77) 
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1 / 1 \ 2 k 2 

— I — - 1 (ki • k 2 )(ki2 • k 4 )— - G bb (k u k 2 , fc 12 , fc 4 ) + 23 perm. . 
2 5 \cj J fci2 



(7? 



1 / 1 \ 2 A; 2 

^5 f — - 1 J (ki2 ■ k 2 )(k 3 ■ k 4 )T J - G bb (-kv2, k 2 , fc 3 , fc 4 ) + 23 perm. . 



'12 



(79) 



1 (" 

2 4 Vc? 



2 7.2 



k 2 k 
^j ( 4r - 1 ) (ki2 • k 2 )(ki2 • k 4 )-^ 3 - 1 G bb (-k 12 , k 2 , k 12 , h) + 23 perm. . 

fc 12 



where 

kia = kx + k 2 , M = k 3 + ki + k 12 , AT = fci + k 2 + fc 3 + fc 4 
and F, G ab , G ba , G bb are defined as follows: 

F(a 1 ,a 2 ,m) 

= — - r2aia 2 + («i + a 2 )m + "^ 2 1 , 

m J L J 

G a6 (ai,a 2 ) 

= jt^ [2a x a 2 + (ai + a 2 )M + M 2 ] 

3 r / xn 12 

,,n^„ [2aia 2 + (ai + a 2 )M\ + otia 2 , 



M 2 ^ 4 



MK 5 



G f 6a(ai,a 2 ) 

" A^K + M^ (ai + « 2 + M ) + W]CS [ 2 «i«2 + 2( ai + a 2 )M + M 2 } 
3 12 

+ TTTT77 [2«1«2 + («1 + ot 2 )M] + ,„^ ttia 2 



M 2 ^ 4 



MK 5 



G ? 66 (ai,a2,a 3 ,a 4 ) 

[2a 3 a 4 + (a 3 + a 4 )M + M 2 ] 



+ 



+ 



M 3 K 

1 

2 
M 3 K 3 

6 

M 2 K 4 



2a 3 a 4 («i + a 2 ) + (2a 3 a 4 + (ai + a 2 )(a 3 + a 4 )) M + ^ «jM 

i=i 

4 

2 JJ a« + (2a 3 a 4 («i + a 2 ) + aia 2 (a 3 + a 4 )) M + ^J ajOjjM 2 

i=l i<j 

iW)(^£i) + -j±,n 

j=l / \ i=l V i=l 



a,- . 



iO) 



(82) 



53) 



i4) 



55) 



Note that in G ab , G ba and G^, the K and M are defined as K = k\ + k 2 + fc 3 + fc 4 and 
M = fc 3 + ki + &i2 (which changes correspondingly in permutations), but not in terms of 
a,-'s. 
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